Interference Alignment with Quantized 
Grassmannian Feedback in the K-user MIMO 

Interference Channel 

Mohsen Rezaee, Maxime Guillaud 
Institute of Telecommunications, Vienna University of Technology, 
Vienna, Austria - e-mail: {guillaud, mohsen . rezaee}@tuwien . ac . at 



Abstract — A simple limited feedback scheme is proposed for 
interference alignment on the K-user Multiple-Input-Multiple- 
Output Interference Channel (MIMO-IC). The scaling of the 
number of feedback bits with the transmit power required 
to preserve the multiplexing gain that can be achieved under 
perfect channel state information (CSI) is derived. This result is 
obtained through a reformulation of the interference alignment 
problem in order to exploit the benefits of quantization on the 
Grassmann manifold, which is well investigated in the single-user 
MIMO channel. Furthermore, through simulations we show that 
the proposed scheme outperforms the naive feedback scheme 
consisting in independently quantizing the channel matrices, in 
the sense that it yields a better sum rate performance for the 
same number of feedback bits. 

I. Introduction 

Multiple-antenna ttansceivers improve the performance of 
a wireless communication link compared to single-antenna 
systems. The increasing demand for high throughput and 
reliable transmission necessitates efficient use of Multiple- 
Input-Multiple-Output (MIMO) systems. The availability of 
channel state information (CSI) at the transmitter is crucial in 
order to fully exploit the performance improvement of MIMO 
systems. 

In scenarios where the channel is not reciprocal (such 
as frequency-division duplex systems) the CSI has to be 
quantized and fed back to the transmitter. The mismatch 
between the true channel and the quantized channel results in 
a degradation in performance. This degradation is more severe 
when it comes to multi-user systems because the channel 
mismatch not only reduces the effective channel gain but also 
causes interference for the other users. 

Interference alignment (IA) is a precoding method that 
achieves the optimal multiplexing gain (also called the degrees 
of freedom, DoF) over the if -user interference channel when 
perfect CSI is available at the transmitters [1J. This method 
designs the precoders such that the total interference at each 
receiver lies in a space with minimum dimensions so that 
the rest of the dimensions can be used for interference-free 
decoding. I A has been introduced for the A'-user MIMO IC 
in 0. 

Extensive research has been made on limited feedback 
schemes for single-user MIMO systems J5] and references 
therein]. In J4), codebook design is investigated when the 
receiver selects the best unitary precoder from a finite code- 
book and feeds back the index of the selected precoder to 



the transmitter. [4| shows that the optimal design for such a 
codebook is equivalent to the Grassmannian subspace packing 
problem. Some useful quantization bounds on the Grassmann 
manifold are derived in 0, J6). In Q, quantization of the 
precoding matrix using a random vector quantization (RVQ) 
codebooks is investigated which provides insights on the 
asymptotic optimality of RVQ codebooks. 

In the interference channel the problem is explored in the 
context of interference alignment over frequency selective 
channels for single-antenna users in |8| and for multiple- 
antenna users in [9|. Both references provide DoF-achieving 
quantization schemes and study the required scaling of the 
number of feedback bits. From another point of view, [10] 
provides an analysis of the effect of imperfect CSI on the 
mutual information of the interference alignment scheme. The 
authors in Hill proposed a method to reduce the quantization 
error w.r.t. the classical scheme; the method involves a 
computationally heavy iterative algorithm which must be ran 
for each codeword and for each channel realization. 

In this paper we present a simple feedback scheme 
based on a Grassmannian representation of CSI and we 
characterize the proper scaling of the number of feedback bits 
with the transmit power in order to preserve the multiplexing 
gain achievable under perfect CSI. Using the structure 
of the interference alignment equations, we remove the 
redundant information in the channel quantization procedure. 
Simulations are also provided to compare our proposed 
method to existing feedback methods. 

The remainder of the paper is organized as follows. In 
Section HI] the system model is described. The reformulation 
of the interference alignment problem is provided in Section 
Hill The limited feedback method is presented in Section |IV] 
The achievable rate is analyzed in Section [V] Simulation 
results are presented in Section |VTJ and conclusions are drawn 
in Section IVHI 

Notation: Non-bold letters represent scalar quantities, bold- 
face lowercase and uppercase letters indicate vectors and 
matrices, respectively. Ijv is the N x N identity matrix. The 
trace, conjugate, transpose, Hermitian ttanspose of a matrix or 
vector are denoted by tr(-), (•)*, (-) T , (-) H respectively. The 
expectation operator is represented by E(-). The determinant 
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of a matrix (or absolute value of a scalar) is represented by 
| • |. The Frobenius norm of a matrix is denoted by || • ||p 
while the two-norm (spectral norm) of a matrix is represented 
by || • H2. A block diagonal matrix is denoted by Bdiag(-) 
with the argument blocks on its diagonal. CAf(0, 1) denotes 
the complex Gaussian circularly symmetric distribution with 
zero mean and unit variance. Finally, X max (-) denotes the 
maximum eigenvalue of a matrix and log represents the logari 
thm in base 2. 

II. System Model 

An interference channel is considered in which K transmit- 
ters communicate with their respective receivers over a shared 
medium. For the sake of notational simplicity, we consider 
the symmetric case where each transmitter has M antennas 
while each receiver is equipped with N antennas, although 
the method applies to more general non-symmetric settings as 
well. 

Each transmitter employs a precoding matrix to transmit d 
data streams to its respective receiver. The vector at receiver 
i reads 

yi = H,. V.x. + ^ 11 . V x . + »< (1) 

l<j<K 

in which Hy € ^NxM j s me cj^jjjjgi ma trix between 
transmitter j and receiver i, "Vj £ c Mxrf anc j Xj . g ^y 1 are the 
precoding matrix and the data vector of transmitter j, respec- 
tively. Furthermore, ~ CAf(0, <t 2 In) is the additive noise 
at receiver i. Assuming E [xjX?] = -jl<j, j = 1, • . • , K and 
using truncated unitary precoders, the transmit power for each 
user is equal to P. The transmit signal-to-noise ratio (SNR) 
is defined as 

III. Proposed Grassmannian Feedback Scheme for 
Interference Alignment 

Let us consider the interference alignment problem of 0, 
in which the channel state information is fed back from the 
receivers to the transmitters and the transmitters find the 
precoding and projecting matrices that satisfy the alignment 
conditions based on the received feedback. This scenario arises 
when the receivers cannot share their channel estimates. The 
ith receiver estimates the channel matrices Hy, j 7^ i and 
feeds back the necessary information to all the transmitters 
so that every transmitter is capable of solving the alignment 
problem. Alternatively, one could consider a single central 
processing unit where the IA solution could be computed and 
subsequently distributed to the transmitters. This distinction is 
immaterial, as the results presented here apply to both cases. 
In this section we consider perfect CSI feedback in order to 
to highlight the intuition behind our limited feedback scheme. 

Assuming that global CSI is available at a given location, 
the precoders Vj, i = 1 ... K can be designed to align the 
interference at each receiver into a N — d dimensional space 
Therefore, interference-free transmission for d streams per 
user becomes possible for a feasible setting. We will further 
assume that (K — 1)M > N, which represents the cases 
of interest where interference would occupy all dimensions 



of the receive subspace in the absence of alignment. A 
solution to the IA problem exists (see [12] and more recently 
lPT3l . lfl4l for feasibility criteria - here we will assume that 
the dimensions and the considered channel realizations are 
such that the problem is feasible) iff there exist full rank 
precoding matrices Vj, j — 1,...,K and projection matrices 
Uj £ C Nxd , i = 1, K such that 

l I'H V. Vi,j£{l,...,K}, j^i, 

rank(Uf H^Vi) = d. 

The first condition can be written from the point of view of 
receiver i in the form 

U?H<V_< = Vi £ {1, ...,#} (3) 

in which V_; = Bdiag(Vi, . . . , Vj_i, V i+ i, . . . , V K ) and 
Hi = [Ha, ...,Hii_i,Hi i+ i, ...,HiK] which is a N x (K - 
1)M matrix. We consider the thin QR decomposition of Hf, 
defined as Uf = F l C l in which Fj is a (K - 1)M X N 
truncated unitary matrix and Cj is a N x N matrix. Note that 
in many cases of interest, such as when the channel matrices 
have coefficients independently drawn from a continuous dis- 
tribution, Hi is full row rank with probability 1, and therefore 
C, is invertible almost surely (a.s.). If we employ the filter 
(C" 1 ) 11 at the ith receiver before projecting the signal into 
the interference-free subspace then we can solve the alignment 
problem at the transmitters based on Ff instead of H;, since 

Uf Ff V_, = Uf Cr H cf if V_ 4 = 
=> (Cr 1 U. i ) H H 4 V_. i = 0. 

This indicates that if the transmitters design 

({Vi}iLv{Ui}E=i) for { F F}£i' tnen we can suppress the 
interference at receiver i using the Vi's as precoders and 
(C" 1 !^) 11 as the projection filter at receiver i. 

The importance of this observation is that F^ is a truncated 
unitary matrix and its structure offers a reduction in the amount 
of required feedback. We now show that the knowledge of the 
column space of F; is sufficient to compute precoding matrices 
that fulfill the alignment conditions. 

Considering perfect CSI feedback on the Grassmann manifold, 
the feedback can be considered to take the form of a matrix 
whose columns span the same subspace as those of Fj, and 
which can be written as F^Ui with U,; being a unitary matrix. 
The transmitter finds the alignment solution based on FjUj. 
Therefore we have 

Uf (F,U l ) H V_. i = (U 4 t%) H Ff V_ 4 = 0. (5) 

Letting Uj = U;Ui and using the same argument as in 
dU, the interference can be completely removed by using 
the filter (C^U,) 11 = (C^LW,) 11 at receiver i. This 
means that the receiver should rotate the alignment projection 
matrix (Uj) using the same rotation employed for feedback 
and then transform it by C^ 1 . Assuming that the process 
of reconstructing F;Uj from the fed back message is 
deterministic and known to the receiver, F^Ui and therefore 
Uj can be assumed to be known at the receiver. Uj, however, 
is not known at the receiver and must be communicated or 
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evaluated through other means. 

It is obvious from © that the CSI feedback problem in this 
context is equivalent to feeding back a point (represented by 
the truncated unitary matrix Fj) on the Grassmann manifold 



reads 



*(K-l)M,N 



for each of the K users. In the next section, we 



consider the quantization of this CSI information. 



IV. Quantized CSI Feedback Scheme 

In this section we consider the case where the alignment 
equations are solved based on the (error-free) feedback of 
a quantized version of the CSI, using the Grassmannian 
representation outlined in the previous section. Let us assume 
that receiver i perfectly knows its channels from all interfering 
transmitters and therefore can perform the QR decomposition 
Hf = FiC,;. Receiver i quantizes the subspace spanned 
by the columns of F^ using Nf bits and broadcasts the 
index of the quantized codeword to the transmitters or any 
other processing unit in charge of computing the V,'s. We 
further assume that the transmitters and receivers share a 
predefined codebook S — {Si, S 2 « / } which is composed 
of 2 t truncated unitary matrices of size (K — 1)M x N 
and ideally is designed using Grassmannian subspace packing. 
The quantized codeword is the closest point on the Grassmann 
manifold, i.e., 



arg min d c (S x ,Fi 



XX 



H 



YY 



H I 



in which d c (X, Y) = -t= 
distance between X and Y S G(k-i)m,n [ 1 5 1 . 



(6) 

is the chordal 



The interference alignment problem is then solved based on 
{F?}f =1 to find ({V^JU,}^) fulfilling 



UfFfV_ 4 = 0, Vie {!,..., K}. 



(7) 



At receiver i, let us define the (not necessarily unitary) 
matrix = FfP,;. Furthermore, inspired by the perfect 
feedback situation, we define the total receive filter as Gf = 
(C^UiUi) 11 = (C^ 1 Fj i FiU i ) H . It should be noted that 
when Fi lies in the subspace of one of the codebook elements 
(i.e., d c (Fi,Fi) — 0) then JJi becomes a unitary matrix and 
will be equal to U, from Section III. 

Here we show that our limited feedback scheme preserves 
the total multiplexing gain of the channel provided that the 
number of feedback bits grows with the transmit power, and 
we characterize the required growth rate. 
After applying the receive filter Gf to ([TJ, the interference 
leakage (due to imperfect CSI) at receiver i can be written as 



l<j<K 



(8) 



Therefore the power of the interference leakage at receiver i 



U = tr(E x (e 4 ef )) 



K 



p 



(9) 



tr(GfH 4 V_. i V H , t HfG l 



Substituting Gf = (Cr 1 Ff F.U,) 11 and H, = Cf Ff gives 

(10) 



^ = l r ||U?F«F i Cr H C < H FfV_ i ||| 



-r 

^lUfFfF.FfV^Hl. 



Using the alignment equation in (|7]l and the fact that Ff Fj = 
Ijv, yields UfFfFjFfV-i = 0, therefore ([TOj can be 



rewritten as 



Li = -||UfFf (F<Ff - PiF? )V_i|||. (11) 



Using the facts that ||X|| F < v /rank(X).||X|| 2 , ||X|| 2 < 
||X|| F and ||XY|| 2 < ||X|| 2 ||Y|| 2 , we have 



Li = |||UfFf(F 4 Ff - F i Ff)V. 
< P||U?Ff (Filf -FiBf)V_i| 



1 1 2 

i I If 



^PllUflllllFflllKFiFf-F^lllllV-il^ (12) 
= P||(F i Ff-F i Ff)||^ 
<P||(F<Ff -FiFf)\\ 2 F 
= 2Pd 2 c (F l ,F t ). 

The second equality holds because Uf , Ff and each block 
of V_i are truncated unitary matrices which implies that the 
two-norm of these matrices is 1. 

Using [5, Corollary 1] and lfT31 Theorem 5], if the codebook 
is generated using the sphere-packing procedure, the maximum 
value of quantization error in terms of the chordal distance can 
be upper bounded as 

2 

(13) 



(02^)^0 

in which A imax = max F , e g (if _ 1)M N 4(F i; Fj), the constant 
c is the coefficient of the ball volume in the Grassmann 
manifold which is given in ||5] Theorem 1] and N& = 
2N((K — 1)M — N) is the real dimension of G(k-~l)m,n- 
Therefore from ( fl2] >. (Tl3T > the leakage term can be upper 
bounded as 



Li < 2PA4 ax < 



8P 



(c2 N f) 



(14) 



In order to achieve the same multiplexing gain as in the 
perfect CSI situation we have to bound the interference leakage 
with a constant value independent of the transmit power while 
letting P — > oo. From ( fl4l i. it is obvious that this requires that 



4 



2 n g scales at least linearly with P. Therefore, we can bound 
the power of interference leakage by letting Nf grow as 

Nr 

N f = -±- logP = N((K - l)M - N) log P. (15) 

Therefore, at receiver i, after applying the receive filter Gf 
to the received signal (HJ we have 



*y% HjjXj 



n, 



in which y z = Gfy l , H l4 = GfH ll V l and n 4 = Gfm + e; 
which has a bounded power independent of P. Note that 
inter-stream interference is not considered here, since we 
assume that Hjj is perfectly known at receiver i. 

It should be noted that the receive filter Gf can induce a 
reduction in the power of the signal of interest. However this 
will not affect the multiplexing gain, since it can be shown 
that 1 1 1 1 1. is proportional to ||Ff F. 
be bounded below as 



F and this term can 



|F?Fj 



\ Z F =N-di(Fi,Fi) >N- A 
and the desired power is written as 

E(tr(H 4l x,xfflH)) = P||H 



2 



(17) 



(18) 



Since A imax scales with P 



the reduction in the desired 
signal power will be a constant value independent of P. 
Therefore we can conclude that the desired signal power is 
of the order of P. Furthermore, it is clear from the definition 
of the d x d matrix Hjj that it is almost surely full rank, since 
Hjj is independent of Gf and Vj. Therefore, we conclude 
that our scheme can achieve d DoF per user if the quantized 
CSI scales according to dlSi . 

V. Achievable Rate Analysis 

In this section we analyze the sum rate achieved by the 
proposed quantization scheme. In particular we bound the rate 
loss incurred by the quantization to find a lower bound for 
the achievable sum rate. Here also we will observe that the 
rate loss would be a constant value independent of the power 
provided that the number of feedback bits is scaled according 
to (fT~5T > which reaffirms the achievability of the total DoF. 

A. Rate Loss Due to CSI Quantization 

Let us first consider the rate P* achievable by user i (after 
filtering by Gf at the receiver) with finite rate feedback: 



K = log 



Qi: 







- log 


id + 4qi 




(T 



(19) 



in which Q| = Gf H u V t Vf H^G, and Qj = 
Gf HjV.jV^jHf Gj. For perfect CSI, the interference is 
totally removed and ( fT9l simplifies to 



r; = log 



(20) 



The rate loss due to quantization reads 



AR l = R l p -R\ 



= log 
< log 



i d + ^Ql 
Id + ^Qi 



.4 



(21) 



(16) where A = log |l d + £(Q l s + Qj)| - log |l d + ^Q|| and 



the inequality follows from the fact that QJ is positive semi- 
definite and therefore A > 0. 
AR 1 can be further bounded as 



AR l < log 



^2 ^1 



<dlog(A max (I d + ^Qj)) 
rflog(l + -^ Amax (Qj)) 



(22) 



dlog(l 



P 



a 



where we have substituted for the interference covariance 
matrix as PQ\ = PGPH,:V_iV H ,H?G,: and used the 



fact that ||X 
HG?H,;V_j| 



= A max (XX H ). From eqs.(l9b— (fT2l> we have 
< 2^(Fj,Fj) which gives 



2P 

AP* < dlog(l H 2"6^(Fj,Fj)). 



(23) 



From d23l it is clear that the rate loss would be constant if 
the number of feedback bits is scaled according to ( [13] ) which 
indicates that full DoF is preserved. 

Another observation is that when the number of feedback 
bits is scaled according to dl~5b . the upper bound of the rate loss 
in (fJTJ will not change by increasing the power however the 
bound gets tighter as power increases. The difference between 
the upper bound and the actual rate loss in ( f2TI ) when P — > oo 
can be written as 



lim A= lim log |Lj + PQ\(a 2 I d + PQs) _1 | 



lim log 

P->oo 



(24) 



which goes to zero since Qg is full rank almost surely and 
when feedback scales according to ( fl~5l ) we have ||Q|||2 — > 
for P — > oo. Therefore the bound in (1211 gets tight. This will 
help us to get a closer bound on the sum rate achievable by 
the quantization scheme as we will see later. 

Using the bound derived for the rate loss in d23l . we can 
lower bound the rate achievable by the quantization scheme 
as follows 

pi = pi - ap 4 



op (25) 
rflog(l + — ^(Fj,Fj)). 



Note that the bound in ((25) holds for any particular codebook 
and also for any realization of the channel matrices. We will 
now consider the average loss in the achievable rate when the 
channel and/or the quantization codebook are random. 
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B. Average Rate Loss under Random Vector Quantization 

If RVQ is used, the relevant performance metric is the 
average sum rate over all possible codebooks and channel 
realizations. Averaging the chordal distance over the channel 
matrices is equivalently performed by averaging over Fj. 
Therefore the average rate for user i reads 

op 

Es,h(J$ > E H (i%) -dE 5)F< (log(l + — ^(F^F,))) 

2P 

> Eh(J$) - dlog(l + E SiF< (^(Fi, Fi))) 

(26) 

where the second inequality follows by application of Jensen's 
inequality to the log function. The term Es ) p 1 (d^(Fi, F^)) 
represents the expected value of the distortion while using 
a random codebook, and can be further bounded using |15 
Theorem 6], which can be summarized as follows: for asymp- 
totically large codebook size, when using a random codebook 
for quantizing a matrix F arbitrarily distributed over the 
Grassmann manifold, the fc-th moment of the chordal distance 
L)( fe ) = E S)F (d^(F,F)) to the quantized point F can be 
bounded as 



(N G + k)(c2 N f)~^ 



< D (k) < 



(27) 



where T(-) denotes the Gamma function. Using the upper 
bound in (l27l i for k = 2, the average rate of user i can be 
lower bounded as 



Es,H(iC) > E H (J?1) - dlog 1 



IP 



r(£) 



^(c2 w /) : 



(28) 



C. Per-User DoF in a Multi-User System 

An interesting consequence of the rate-loss analysis con- 
ducted previously can be observed when each receiver uses 
its own scaling of the CSI quantization codebook size with P. 
Formally, let Nl denote the number of bits used by receiver 
i to quantize F^. Furthermore, assume that Nj scales with P 
such that 



lim 



N 



f 



(29) 



p->oo N((K — 1)M — N) log P 
exists and is finite. Then DoF achievable by user i is given by 



d% = d% min (pa, 1) 



(30) 



where is the achievable DoF of this user with perfect CSI. 
The proof follows simply from ( 1251 ) by taking the limit of 
the lower bound when P — > oo, and therefore holds both for 
RVQ and for codebooks generated using the sphere-packing 
procedure. 

Practically, this means that the DoF achieved by each user 
is independent of the quality of the feedback provided by the 
other users. This result, obtained here for precoding based 
on the IA criterion, is consistent with the scaling obtained in 
1 16 1 for centralized schemes using different precoding schemes 
such as zero-forcing. 



VI. Simulation Results 

A. Performance results using RVQ 

In this section, the performance of the proposed scheme 
is evaluated through numerical simulations. The performance 
metric is the sum rate evaluated through Monte-Carlo simu- 
lations employing the precoders designed using interference 
alignment. 

The achievable sum rate of the MIMO IC using interference 
alignment precoders under the assumption that the input 
signals are Gaussian can be written as 



K 

Rsum = E l0 £ 



K 



a 2 I N + £HyV J V?Hg 



K 



K 



(31) 



A three-user IC with two antennas per node and one data 
stream for each transmitter is considered. Entries of the 
channel matrices are generated according to CJ\f(0, 1) and 
the performance results are averaged over the channel realiza- 
tions. Our proposed method is compared to the naive method 
where the interfering channel matrices toward a receiver are 
independently vectorized and quantized based on the idea of 
composite Grassmann manifold J9) and finally the indices of 
the quantized vectors are fed back to the transmitters (denoted 
by "Naive" in the figures). 

Figure Q] shows the achievable sum rate versus transmit 
SNR for Nf = 5 and 10 feedback bits when the precoders 
are designed based on the quantized feedback. Instead of the 
optimal subspace packing codebook, a random codebook is 
used where in our proposed method the codebook entries 
are independent (K — 1)M x N random truncated unitary 
matrices generated from the Haar distribution. For the naive 
method, random unit norm vectors are used in the codebook 
construction. Clearly the proposed scheme outperforms the 
naive method for the same number of feedback bits. 
The sum rate in OTT l is achievable when optimum receivers 
(not including the projection filters Gf) are used at the re- 
ceivers. Since the achievable scheme in Section|IV]is using the 
projection filters Gf, we evaluated the performance achieved 
by this scheme, defined as 



R' 



K 



sum / t 

i=l 

K 

E l0 § 



o»G?G t 



K 

log cr 2 Gf Gj + ^2 Gf H,- 

K 

Gfll; , V, V '' Hj- ' C; , . (32) 
Results are provided in Figure [2] 

The slope of the curves on Figure [2] at high SNR gives us 
an indication of the DoF. It is clear from Figure [2] that the 
slope of the rate curve with quantized feedback matches that 
of perfect CSI when the number of feedback bits is scaled 

according to (B). We have used Nf = [0, 7, 13, 20, 26] bits 

2N f 

and their corresponding powers (P = 2 n g ) to generate the 
curve which exhibits the achievable DoF. Simulations were 
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Fig. 1. Sum-rate comparison of quantization methods, for the 
3-user MIMO IC, N = M = 2 using OTJ. 
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Fig. 2. Sum-rate according to (1321 of the proposed method for 
different number of bits, for the 3-user MIMO IC, N = M = 
2. 



performed only up to SNR of 20 dB due to the complexity 
associated to the exponential size of the codebook. 

B. Perturbations on the Grassmann manifold 

In order to demonstrate the achievable sum rate according 
to ( TTBT l for high power region (which corresponds to high 
number of feedback bits) we replace the quantization process 
with a perturbation which approximates the quantization error. 
As will be seen, this approach gives us a good approximation 
of the sum rate specially when the number of feedback 
bits increases, while avoiding the prohibitively complex 
simulation of RVQ schemes for very large codebooks. 

Consider a n x p truncated unitary matrix F, representing 
a point on Q n . p . Clearly we can assume p < ^ (otherwise, 
the same point can equivalently be represented by a matrix 
spanning the null space of F, which is of dimension lower 



than S). According to (6), any other point on Q n p can be 
represented in the basis constituted by the columns of a unitary 
matrix W having F as its first columns, i.e. W = [F F c ], as 



F = W 



C 
S 

n—2p 



(33) 



for some F c spanning the null space of F and 



C = 



COS fi 







, S = 



sin Oi 



up _ 

(34) 



for some 6\, 9 p . In this basis, F is trivially obtained as 



W 



I, 



n—p 



(35) 



The squared chordal distance between the two points is 



= d 2 c (F,F) =^>m 2 6, 



(36) 



In order to generate random perturbations of a certain chordal 
distance y/r^ from F one can generate random values for the 
angles 6x,...,9 p such that X)f=i sm2 @i = r o an d compute the 
perturbed matrix from d33l ). for some randomly chosen F c . 

The histogram of the squared chordal distance between the 
actual Fj matrices and their quantized versions Fj obtained 
through simulation of the RVQ quantizer suggests that the 
Gaussian distribution is a good approximation for the prob- 
ability density function. The parameters of this distribution 
can be obtained from the moments in d27b . Note that d27l ) 
only provides bounds on D^ k \ however since those bounds 
are asymptotically tight when the codebook size increases, we 
arbitrarily choose to use the upper bound as an approximation 
of £>(*), i.e. 



rc-M 



No 
2 



D {2) 



(c2 N f) 



is the average and 



D (4) _ p(2) )2 



(37) 



(38) 



is the variance. We propose generate the values for r 
according to the this Gaussian distribution truncated to 
the interval (0, A^ lax ). The process of generating random 
perturbations is summarized in Algorithm 1. 

The sum rate performance of the perturbation method is 
plotted against the actual quantization scheme in Figure [3] 
We use the upper bound in ( f27T > to evaluate the parameters of 
the approximated distribution and the difference in the results 
is negligible when we use the lower bound. It is clear that this 
perturbation method effectively approximates the quantization 
process when the desired performance metric is the sum rate 
specially at regimes that interference is not dominant. This 
allows us to rely on the curves resulting from this method 
to analyze the sum rate for high power regime (in particular 
when the number of feedback bits are scaled according to 
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Algorithm 1 Generating random perturbations around F 

• Draw a random realization of the chordal distance r from 

A^(r-a r 2 ) 

• If r < or r > £^ max then generate a new sample 

« Generate a random vector s = [s±, ■ ■ ■ ,s p ] such that < 
s l < 1 and J2i=i Si = r 

• Compute the angles as 0, = sin -1 ^fsl and construct C 
and S according to ( l34l 

• Generate a random orthonormal basis F c for the null 
space of F 

C 

s 



Compute F = [F F c ] 




u 

I 

s 

3 

(75 



20 

SNR [dB] 

Fig. 3. Comparison of the perturbation scheme from Sec- 
tion IVI-BI (solid) to the real quantization method (dashed), 
for the 3-user MIMO IC, N = M = 2. 



( fl~5T > and RVQ has impractically high complexity) to observe 
the achievable DoF. 

Figure |4] shows the sum rate performance using the pertur- 
bation method compared to perfect CSI and the lower bound 
derived in (1281 . The slope of the sum rate at high SNR regime 
is identical to that of perfect CSI for the perturbation method 
(which accurately approximates the actual sum rate) and also 
for the lower bound derived in (|28T >. As expected earlier, we 
observe that the lower bound gets closer to the actual sum 
rate when the interference power is bounded and the desired 
power increases and with a constant number of feedback bits 
the lower bound becomes looser as the power increases. 

VII. Conclusion 

A new CSI feedback scheme for interference alignment on 
the K-user MIMO interference channel was proposed consist- 
ing in a parsimonious representation based on the Grassmann 
manifold. We characterized the proper scaling of number of 
feedback bits with transmit power in order to preserve the 
multiplexing gain achievable using perfect CSI. Simulations 
results confirm that our scheme provides a better sum rate 



N 

X 



3 

CO 



40 



30 



20 



10 



N f 



25, lower bound 
25 perturbation 



■ log P, lower bound 

' Nf = —If- log P perturbation 
- Perfect CSI 




20 

SNR [dB] 

Fig. 4. Sum rate performance using the perturbation method 
compared to perfect CSI and the lower bound derived in d28l ), 
for the 3-user MIMO IC, N = M = 2. 



performance compared to the naive feedback scheme for the 
same number of feedback bits. 
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